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Abstract 


This paper addresses stabilization and control issues in autonomous capture and 
manipulation of non-cooperative space objects such as asteroids, space debris, and 
orbital spacecraft in need of servicing. Such objects are characterized by unknown 
mass-inertia properties, unknown rotational motion, and irregular shapes, which 
makes it a challenging control problem. The problem is further compounded by 
the presence of inherent nonlinearities, significant elastic modes with low damping, 
and parameter uncertainties in the spacecraft. Robust dissipativity-based control 
laws are presented and are shown to provide global asymptotic stability in spite of 
model uncertainties and nonlinearities. It is shown that robust stabilization can be 
accomplished via model-independent dissipativity-based controllers using thrusters 
alone, while stabilization with attitude and position control can be accomplished 
using thrusters and torque actuators. 


1 Introduction 

The proposed NASA asteroid redirect mission has created considerable excitement 
and interest in the public as well as in the worldwide science and engineering com- 
munities. A study by the Keck Institute [1] has concluded that it is feasible to 
autonomously capture and return an entire 7-m diameter, 500,000-kg near-Earth 
asteroid to a high lunar orbit. An alternate approach, consisting of picking a boulder 
off of a larger asteroid, is also being studied. Investigations are in progress at NASA 
and other organizations to conceive and develop alternative approaches and methods 
for capture, manipulation, and transport of asteroids. 

Another important challenge that is technically similar to asteroid capture is 
orbital debris mitigation. The near-Earth as well as geostationary orbital debris 
population is continuously growing and its growth is expected to continue in the 
future due to ongoing space activities. On-orbit satellite explosions and collisions 
(accidental or intentional) create even larger numbers of debris items. Space debris 
poses a serious threat both to human-occupied vehicles and to commercial satellites. 
Some suggested approaches for orbital debris mitigation would involve docking a 
capture spacecraft with the debris or making a physical impact. This is a techno- 
logically challenging task as most debris are non-cooperative and possibly tumbling 
with unknown spin, precession, nutation, amplitude changes etc. 

A third technically similar challenge is autonomous on-orbit capture and servic- 
ing of defunct or functioning satellites. Eor example, on-orbit refuelling or repairs 
can extend a satellite’s life at a fraction of the cost of constructing and launching a 
new satellite. Also, with the rapid ongoing advancement of micro-component tech- 
nologies, it may be desirable to replace components of older functioning satellites 
with advanced components and significantly increase the capabilities. (Eor refuelling 
or refurbishment of functioning satellites, the capture and manipulation should be 
somewhat simpler since they would likely be cooperative objects). 

A related problem is asteroid strike threat mitigation, which would involve chang- 
ing an asteroid’s trajectory to avoid collision with Earth. 
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The problems described above, i.e., asteroid capture, space debris mitigation, on- 
orbit spacecraft servicing, and some methods for asteroid strike threat mitigation, 
involve capture and manipulation of non-cooperative space objects (NCSOs) and 
would require a common set of technologies. 

The main technical problem in capture and manipulation of NCSOs is that they 
have unknown mass-inertia properties and their dynamic state (such as angular ve- 
locity vector- spin, precession, nutation) are unknown. Additionally, asteroids have 
irregular shapes and composition (e.g., ranging from solid monolithic bodies to col- 
lections of loose debris) . Methods for capture and manipulation of non-cooperative 
space objects would likely require capture spacecraft having large solar arrays, multi- 
link manipulators, deployable bags, and/or other capture mechanisms. In addition to 
uncertainties in the non-cooperative space object’s properties and state, the capture 
spacecraft will likely have significant lightly damped elastic modes with uncertain 
and varying modal parameters (frequencies and mode-shapes). Because of these 
problems, capture, stabilization, and post-capture control of non-cooperative space 
objects constitute a formidable technical challenge. 

This paper addresses some basic issues in capturing a non-cooperative space ob- 
ject and stabilizing the combined spacecraft /object body during and after capture. 
The capture spacecraft is assumed to be a generic multi-link manipulator having 
a branched geometry configuration that includes a relatively large central body to 
which multiple articulated branches are attached through single degree-of-freedom 
(DOF) revolute joints. It is assumed that the object’s spin/tumble rate is relatively 
small, so that pre-capture co-spin of the spacecraft (i.e., matching the object’s angu- 
lar velocity vector) is not performed. It is also assumed that orbital dynamics do not 
significantly impact the problem addressed in this paper, which will be the case when 
the time duration of the capture maneuver is small compared to the orbital period. 
For applications where orbital dynamics can significantly impact the system motion 
(e.g., in low Earth orbits) the orbital dynamics will also have to be considered. This 
paper presents in detail and extends the preliminary analyses and results presented 
by the authors in [2] for which there was no written version. 

1.1 Control Objectives 

Prior to capture, it will be necessary to position the spacecraft at a desirable posi- 
tion relative to the NCSO, with a desired posture (defined by the base body attitude 
and the link joint angles) so that the end-effectors are appropriately positioned for 
capture. During capture, it will be necessary to stabilize the spacecraft /NCSO com- 
bination in the presence of unknown contact forces acting on the end effector. In 
the post-capture stage, it will be necessary to stabilize the combined body, i.e., to 
bring the rotational and translational motion as well as all elastic motion to rest. 
It is assumed that the primary actuators on the spacecraft would be thrusters, pos- 
sibly augmented by torque actuators such as control-moment-gyros (CMGs). It is 
also assumed that each thruster can produce the commanded variable thrust (i.e., 
proportional thrust, e.g., see [3], [4]) rather than operating in an on/off mode. In 
addition, each link joint is assumed to have a torque actuator. In the pre-capture 
stage, the spacecraft dynamics are known to some extent (although uncertainties 
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will still exist in the mass inertia properties and elastic mode dynamics). In the 
mid-capture and post-capture stages, however, large unknowns and uncertainties 
will exist in all aspects of the combined spacecraft/NCSO dynamics. Hence, ro- 
bust stabilization during and after capture is a major technical challenge. In the 
post-capture phase, the link joints would likely be locked, and the dynamics of the 
combined spacecraft/NCSO can be represented by a single body. During capture, 
the stabilizing control laws (which keep the velocities and elastic motion near zero) 
would be in effect, so that the stabilized spacecraft subjected to contact forces (‘dis- 
turbances’) would likely result in relatively small angular velocities that would be 
eventually brought down to zero after the capture operation is completed. While at- 
titude and posture of the spacecraft (represented by the base body attitude and the 
link joint angles) are important in the pre-capture stage, the primary post-capture 
control objective is to stabilize the combined body, i.e., bring it to rest. A secondary 
objective is to achieve a desired attitude of the combined body in order to prepare 
for transporting it to another location (such as a Lunar Lagrange point in the case 
of asteroid capture). 

1.2 Outline of the Paper 

Section 2 addresses control of single-body spacecraft, such as a post-capture space- 
craft/NCSO combination. Stabilization is addressed first, followed by attitude and 
position control. Both rigid and flexible spacecraft are considered. Section 3 ad- 
dresses stabilization as well as position, attitude, and posture control of multi- 
body /multilink spacecraft. Section 4 presents concluding remarks. 


2 Single-Body Spacecraft Stabilization and Control 

2.1 Stabilization of Rigid Spacecraft 

In the mid- and post-capture stages, the objective is to stabilize the combined space- 
craft/NCSO, i.e., to eliminate the inertial translational and rotational velocities. 
This can be accomplished by using thrusters alone, i.e., without the need for torque 
actuators. The next subsection addresses the case when only thrusters are available. 


2.1.1 Stabilization Using Thrusters 

Suppose there are m bi-directional thrusters at locations ri,V2, ■ ■ ■ Vm, (ri G 1 ?^) 
relative to the center of mass (c.m.) of the combined spacecraft/NCSO (if post- 
capture) , expressed in a body- fixed coordinate system and the ith thruster produces a 
scalar force T) in the direction di where di G is a unit vector. Thus, the ith thruster 
produces the vector force diFi expressed in the body-fixed coordinate system. It 
is assumed that ri,r2,--- ,rm are such that rank[fidi, f2ci2, • • • is 3 , where 

overhead ~ denotes the cross-product matrix of a vector, i.e., if ^ = [^X7^y,^z]'^ 7 


i = 


0 -6 (y 

& 0 
-^y (x 0 


( 1 ) 
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The equations of motion are given by 


M 0 
0 J 


V 

UJ 


+ 


Ai{ujv) 

UjJui 


cl\ C?2 • • • dm 

ridi V2d2 ■ ■ ■ 'l"mdm 


(2) 


where A4 = mo/ 3 , ^0 being the mass (including the NCSO if post-capture), I 3 
denotes the 3x3 identity matrix, J is the 3x3 moment-of- inertia matrix, and u, oj 
are the inertial translational velocity vector of the c.m. and the angular velocity 
vector, expressed in the body-fixed coordinate system; F = [Fi, F 2 , . . . is the 
m X 1 thrust vector. It is assumed that m > 6 because at least 6 bi-directional 
thrusters are needed to simultaneously control rotational and translational motion. 

Denote 


X = M = 


M 0 
0 J ’ 


and r 


d\ d 2 ... dm 
fidi r2d2 ■ ■ ■ f^dm 


( 3 ) 


In the post-capture stage, the c.m. of the combined spacecraft/NCSO is unknown, 
(as are Ai and J), therefore r^s are not known. It is assumed that T continues 
to have full rank (= 6 ) regardless of the unknown c.m. location of the combined 
spacecraft/NCSO. (This can be ensured by analyzing scenarios for a range of NCSO 
mass-inertia properties. It may also be desirable to include sufficient redundancy, 
i.e., deploy more than 6 thrusters as needed). 

Suppose inertial velocity sensors are placed at the thruster locations. (From a 
practical viewpoint, placement of sensors near the thrusters would produce noisy 
measurements and it may be necessary to use noise-reduction filters. However, it 
would not affect closed-loop stability and would likely have a relatively small effect 
on the closed-loop motion). The velocity measurement Vfi in direction di at the 
location of thruster i (ignoring sensor noise) is given by 

Vfi = dj{v + ujri) = dj{v + rJoj) (4) 


Therefore the velocity measurement vector G 3ft™' is given by 

Vr = Vf = T'^X 


( 5 ) 


Consider the control law: 


F = —GrVm where Gr = G^ > 0, 


( 6 ) 


which yields the closed loop equation: 

MX + [{MZjv)'^ {ujJujff + TGrT^X = 0 (7) 

It is shown below that this control law provides closed-loop global asymptotic sta- 
bility. Furthermore, the stability is robust to (0, oo)-sector actuator nonlinearities 
as well as first-order actuator dynamics. 

Theorem 1-The closed- loop system consisting of (2), (5), and ( 6 ) is globally asymp- 
totically stable. 

Proof- Consider the Lyapunov function 

V{X) = X^MX ( 8 ) 
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Because M = M'^ > 0, y is positive definite. Differentiation with respect to (w.r.t.) 
time t and simplification utilizing properties of the cross product yields 

V = -2X'^rGrr'^X (9) 

Because F has full rank, D < 0 (negative definite), i.e., the closed- loop system is 
globally asymptotically stable □ 

Note that, although the system comes to rest (v,uj —)■ 0 as t —)■ oo), the final 
inertial position and orientation (attitude) are arbitrary. An important feature of this 
control law is that the knowledge of the parameters or the c.m. location (mo, J, ri), 
is not needed, and therefore the stability is robust in spite of parameter errors or 
lack of knowledge of the parameters. 

The robust stability property also holds in the presence of (0, oo) sector actuator 
nonlinearities, as well as first-order actuator dynamics, as shown next. 

Suppose the actuators have (0, oo)-sector nonlinearities, i.e., the actual actuator 
force is given by 

F = (f){-Gryr) (10) 

where (p{.) = , 4 > 2 {-) ■, • • • , , and 

a(f>i{a) > 0 for (T / 0 and </>i(0) = 0 for i G [1, . . . , m] (11) 

4>i{.) are assumed to be continuous functions. 

Theorem 2- The control law in (10) provides global asymptotic stability of the 
closed-loop system in the presence of (0, oo)-sector actuator nonlinearities provided 
Gr is a diagonal matrix. 

Proof- Let 

•0j(o-) = -(j)i{-a) for i G [1, 2, . . . , m] (12) 

Then ifi also belong to the (0, oo)-sector. Proceeding with the Lyapunov function in 
Theorem 1 yields 


m 

V = -2X'^rip(Grgr) = -2'^yriGrii'lpi(gri) < 0 (13) 

i 


Thus V is negative semidefinite. Also V = 0 implies yr = 0 (because of (11)), 
therefore X = 0 (because T is full rank). Therefore the closed-loop system is globally 
asymptotically stable. □ 

It can be shown that this control law also retains global asymptotic stability in the 
simultaneous presence of linear time-invariant (LTI) first-order actuator dynamics 
followed by (0, oo)-sector actuator nonlinearities. The proof is along the lines of [5] 
(Theorem 16 in Chapter 2) and is omitted. 

Remark- As stated in Section 1.1, it is assumed throughout this paper that the 
thrusters can produce proportional variable thrusts rather than operating in an 
on/off mode. If the thrusters are operated in an on/off mode (that usually includes 
a deadzone), it is expected that this control law would yield ultimate boundedness 
(e.g., a limit cycle), although further analysis would be needed to investigate that 
case. 
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2.1.2 Stabilization Using Thrusters and Torque Actuators 

When torque actuators are available in addition to thrusters, the equations of motion 
can be written as: 


M 0 
0 J 


V 

U! 


M.OJV 

UjJuj 



xm 



(14) 


wherein it is assumed for simplicity that a single 3-axis torque actuator is used to 
produce control torque Ta G In addition to velocity sensors, suppose angular 
velocity sensors are available. Then the sensor output vector is given by 

yr = [vJco^f = T^X (15) 


where 



Using the control law 

u = —GrVr where Gr = G'^ > 0 


(16) 

(17) 


and following the same procedure as in Sec. 2 . 1 . 1 , it can be shown that the closed- 
loop system is globally asymptotically stable, and that the global asymptotic stability 
is preserved (if Gr is a diagonal matrix) in the presence of ( 0 , oo)-sector actuator 
nonlinearities as well as simultaneous ( 0 , oo)-sector actuator nonlinearities and first- 
order actuator dynamics. 


2.2 Attitude and Position Control of Rigid Spacecraft 

Assuming both thrusters and torque actuators are available, the equations of motion 
using quaternions to represent attitude are; 


MR 
JCo -|- w X Joj 
a 

$ 


rpT TP 

i 


T 

— [cj X a {(3 1) uj] 

_1 T— 

2^ “ 


(18) 


where R denotes the c.m. position in an inertial coordinate system, Tjb denotes 
the transformation matrix from inertial to body-fixed coordinate system, Fi = diFi 
denotes the ith vector thrust expressed in the body-fixed coordinate system, and the 
quaternion parameters satisfy 


a = /3 = 04 — 1; (19) 

a = [«i, a 2 , a^a = l. ( 20 ) 


The control torque r applied to the spacecraft consists of the moments generated 
by the thrusters and the torque Ta provided by the torque actuators such as control 
moment gyros (CMGs). 

r = S™i(ri xFi) + ra (21) 
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Consider the nonlinear attitude control law: 


r = —-^[{a + aiI)Gp + u{l — a 4 )I]a — GrU} 

= ~2 “ 1 “ (/5 + ~ z // 3 /] q : — GrOJ ( 22 ) 

where Gp, Gr are symmetric positive definite (3 x 3) matrices, and v \s & positive 
scalar. The control law for translational motion may be designed using 


Tjslh, h,..., h]F = TfB[di,d2 , . . . , dm]F = (23) 

where are 3x3 positive definite matrices. 

It was shown in [6] that the rotational part of the closed-loop system has exactly 
two equilibrium solutions: [a = a; = 0, /3 = 0] and [a = oj = 0, /3 = —2], which 
correspond to the same equilibrium point in the physical space, i.e., the equilibrium 
state is unique, defined as: R = 0, R = 0, a = 0, cu = 0, and (3 = 0. By following 
the procedure in [6] with a minor modification (i.e., addition of the term {R^ A4R + 
R^G^p°‘'^^R) to the Lyapunov function), it can be shown that this control law globally 
asymptotically stabilizes the physical equilibrium (the origin of the state-space). The 
proof is based on Lyapunov analysis and LaSalle’s invariance principle. 

The required control thrusts and torque can be obtained by using 


■ dl •• 

dm 0sx3 

'F' 

— r 

'F' 



hdi ■ ■ 

fmdm I 3 


— 1 



— ^[(cr {/3 F l)I)Gp — z//3/]q; — Grcv 


(24) 

r is assumed to have full rank, therefore T also has full rank, and the minimum-norm 
solution u is obtained as 


u = 


F 

Ta 


= r'^(rr'^)-^ 


— 2 [(q: + (/? + ^)d)Gp — uj3I]cx — GrUJ 


(25) 


Suppose the inertial position Rp and velocity vp, expressed along the body axes, are 
measured at a point P that is fixed to the spacecraft at a location denoted by the 
vector rp in the body-fixed coordinate system centered at the c.m. Then R and R 
can be obtained as 


R = tJb [R p - rp)] R = Tfp [vp - ujrp] (26) 

for implementation in the right-hand side of (25). 

If there are no torque actuators, r must be produced by thrusters alone, i.e., it 
is necessary to solve the following equation for F 

'T \r^trans d i /^trans r>] 

-p rp _ XL + it\ f 27 l 

~ [-U{^+{P + l)I)Gp-upi]a-Gru\ ^ ’ 

r is assumed to have full rank, therefore the minimum-norm solution F can be 
obtained similar to (25). Note that ri,--- ,rm, and rp (i.e., the c.m. location) 
must be known for computing u in (25) or F from (27). This is possible in the pre- 
capture stage because the spacecraft c.m. location is known reasonably accurately. 
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In the post-capture stage, the c.m. location is not known, therefore it is usually 
not possible to implement this control law. However, while attitude and position 
control are important in the pre-capture phase, the primary objective in the post- 
capture phase is stabilization; therefore attitude and position control may not be 
necessary post capture. If torque actuators are available, they can be used for post- 
capture attitude control (without using the thrusters) for globally stable attitude 
maneuvering using the nonlinear control law (22). 

2.3 Stabilization of Flexible Spacecraft 

The capture spacecraft would likely have significant elastic mode dynamics because 
of long booms, solar arrays, and manipulator links. The problem of stabilization is 
first considered i.e., the translational and angular velocities as well as elastic motion 
of the post-capture spacecraft/NCSO need to be brought to zero. 

2.3.1 Stabilization Using Thrusters 

The equations of motion including elastic motion can be written as: 

At 0 0 i) A4cuv I 3 I 3 ... I 3 

0 J 0 UJ + OjJu = Pi P 2 ... hmxl (28) 

_ 0 0 In, \ [ q \ [ Dq + Aq \ [ _ 

where q{t) G is the modal amplitude; A = diag{tJi, UI 2 , ■ ■ ■ , u>nq) where w* denotes 
the natural frequency of the elastic mode; and D = > 0 {D G i§ the 

damping matrix, pi G represents the location of the ith thruster (including elastic 
displacement) expressed in the body- fixed coordinate system (fixed to the nominal 
rigid body), given by 

Pi = n + ^iq (29) 

where G is the mode shape matrix at thruster i. (This formulation assumes 

relatively small elastic deformations. For large elastic deformations, a fully coupled 
rigid-elastic model would be appropriate). Velocity measurements at thruster lo- 
cations in the thruster direction dj, measured in the local coordinate system (that 
includes local elastic deflections) are: 

Vri = djQi{^iq){v + ujpi + ^iq) (30) 

where 'kj G is the mode-slope matrix at the ith sensor location, '^iq G 

represents the local relative attitude (with respect to the body-fixed coordinate sys- 
tem) due to elastic modes, and 0i(Tig) G is the corresponding transformation 

matrix. Scalar thrusts (denoted by A)) are generated in the local coordinate systems. 
The ith thruster force, expressed as a 3 x 1 vector in the body-fixed system, is given 
by 

Fi = eJ{^^q)diF^ 

8 


(31) 



where .Fj G K is the thrust in the local coordinate system. Equation (28) can then 
be written as 


V Aiujv 

M Cj + ujJuj 
q Dq + Aq 


r n 

h h ••• h F2 

pi P2 ■■■ Pm Q 

L^l ^2 • • • ^mj 

y~ m_ 

L T := BT (32) 


where 

= diag{Q\{'^iq)di, Ql{^2q)d2, ■ ■ ■ ,Q^{'^mq)dmhmxm 

f = [ ^3. • • /3 1 ^ qT (33) 

\_Pl, P2, - ■■ Pm\ 

M = diag(At, J,4J G SRK+6 )xK+ 6)^ @ ^ gfjmxSm^ f ^ J ^ As 

in the rigid spacecraft case, it is assumed that m > 6. It can be shown that the 
velocities at thruster locations in the local coordinate systems are given by 

yr = B'^X = [f^ $]X (34) 

where X = (u^, q^)'^ ■ Note that this system is highly nonlinear and time- varying 

because of cross-product terms and transformation matrices that are functions of 
the modal amplitude. The system is also unknown and uncertain because Ai, J 
are unknown, modal amplitudes q{t) cannot be measured and the modal parame- 
ters (natural frequencies, mode shapes, damping) are never accurately known. In 
addition, the number of modes (infinite in theory) can be very large, resulting in 
unmodeled modes. In spite of these problems, the closed-loop system can be shown 
to be globally asymptotically stable. It is assumed that T has full row rank (= 6) 
for all t. Consider the control law: 

X = -GrVr; Gr = G'^ >0 (35) 

The closed-loop system is given by 

Mujv 

MX+ djJu +BGrB'^X = 0. (36) 

_ Dq + Aq _ 

Theorem 3- The closed- loop system consisting of (32), (34), (35) is globally asymp- 
totically stable, i.e., v{t),Lj{t), and q{t) — )■ 0 as t — )■ 00 . 

Proof- Consider the Lyapunov function 

V{X) = X^MX + q^Aq (37) 

Differentiation w.r.t. t and simplification yields 

V = -2X^BGrB^X - 2q^Dq (38) 
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Thus, V is negative semidefinite. Also, V = 0 only when g = 0 and yr = 0, i.e., v = 
w = 0 (because T has full row rank), T" = 0, and q = 0 (from (32)). Therefore the 
system is globally asymptotically stable. □ 

The significance of this result is that robust stabilization can be accomplished 
using thrusters in spite of time-variation, unmodeled elastic mode dynamics, nonlin- 
earities, and parametric uncertainties. 

Similar to the rigid spacecraft case, it can be shown that the global asymptotic 
stability is preserved in the presence of (0, oo)-sector actuator nonlinearities. How- 
ever, additional analysis will be needed to investigate stability in the simultaneous 
presence of actuator nonlinearities and actuator dynamics. 


2.3.2 Stabilization Using Thrusters and Torque Actuators 

When torque actuators are available in addition to thrusters, (32) is modified as 



Mu)v 



( 03x3 


'F 

MX + 

uJoj 

= 

B 




Dq + Aq 



\ j\ 


'^a 


(39) 


wherein u = [T'^,rJ]^ and it is assumed that three axis torque G 3^^ is produced 
by a torque actuator; G is the mode-slope matrix at the torque actuator 

location, and Qt G is the transformation matrix from body-fixed coordinates 

to local coordinates at the torque actuator location. 

Assuming collocated sensors and actuators, the sensor output pr consists of in- 
ertial velocities uy G at the thruster locations and the angular velocity uit G 
at the torque actuator location, all in local coordinates. It can be shown that 

= B^X (40) 

Using the control law 

U = GrVr (^ 3 .) 

where Gr is positive definite and symmetric, and proceeding as in the previous 
subsection, it can be shown that the closed-loop system is globally asymptotically 
stable, and that this stability property is preserved in the presence of (0, oo)-sector 
actuator nonlinearities when Gr is a diagonal matrix. 



2.4 Attitude and Position Control of Flexible Spacecraft 

When it is required to control attitude and position, it was seen in Section 2.2 
for rigid spacecraft that the knowledge of the thruster locations w.r.t. the c.m. 
was needed. This is generally not possible in flexible spacecraft because of elastic 
motion unless additional instrumentation is used. One approach to position and 
attitude control would be to first achieve the desired translation (c.m. position) 
using thrusters (while ensuring no net moment, i.e., pi x Fi k. 0), and then use 
the torque actuators to achieve the desired attitude. If non-zero moment generated 
due to uncertainties causes significant non-zero angular velocity, it will be necessary 
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to stabilize the spacecraft first (Sections 2.3.1 and 2.3.2) and then perform attitude 
control using torque actuators. Attitude control (using a 3-axis torque actuator and 
collocated attitude and rate sensors) is considered in this section. 

Following [7], the rotational motion can be modeled as 

M{p)p + C{p, p)p + Dp + Kp = B'^Ta (42) 

where p = {uP" , M{p) G 3fj("-9+3)x(nq+3) jg symmetric positive definite ma- 
trix, C{p,p), D G 3fj("-9+3)x(nq+3) corresponds to Coriolis and centrifugal forces and 
damping terms respectively. 

j. _ 0sx3 )l>3xnq ^ _ 0sx3 O 3 X ng 
■^~n ^ ^ ^ ~ r\ n 

yriqXS -^nqXriqj L^n^xS ^nqXUq_ 

From [7], K, D are symmetric positive definite matrices, and 

S {p, p) = M (p) - 2C{p, p) (44) 

is skew-symmetric. The attitude is measured at the torque actuator location and is 
expressed in the quaternion form as in (18). Using the control law as in (22), 

Ta = -^[(a + (/? + l)4)Gp - uj3I]a- GrOJ (45) 

where Gp, Gr are symmetric positive definite and 0 < XuiGp) < 2z/ (where Xm{-) 
denotes the largest eigenvalue), it can be shown (as in [7]) that there are two closed- 
loop equilibrium solutions that are physically identical, i.e., the equilibrium is unique. 
The closed-loop global asymptotic stability can be established by using a procedure 
similar to [7] and the Lyapunov function 

V = p^M{p)p + Kq + GpOi + vfP (46) 

3 Multi-body Spacecraft Stabilization and Control 

In multi-body spacecraft such as space-based manipulators with multiple articulated 
joints, stabilization can be performed by locking the joints first, which essentially 
results in a single-body spacecraft. The methods in Sections 2.1.1, 2.1.2 and 2.3.1, 
2.3.2 can then be used for rigid spacecraft and flexible spacecraft. For position and 
attitude control of multi-body spacecraft, one approach would be to first achieve 
the desired translation (c.m. position) with joints locked, using the thrusters while 
ensuring minimum net moment, and subsequently use the torque actuators and the 
joint torques to achieve the desired attitude and posture. If the translation maneuver 
results in significant angular velocity, stabilization can be performed prior to the 
attitude/posture control. Therefore only attitude/posture control results for flexible 
multi-body spacecraft are presented below. 

Following [7], the flexible multi-body system consists of a central body having 
multiple articulated appendages connected through ne revolute joints. All bodies are 
assumed to be flexible. The dynamics are described by (42) where p = {uP ■, P , 
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0 G K”® is the joint angle vector; M{p) G 3fj(»^9+»^e+3)x(n,+ne+3) jg gynimetric positive 
definite; C{p,p), D G corresponding to Coriolis and centrifugal forces 

and damping terms respectively, are similar to (43). The input vector is 

u = (47) 

where G 3^^ is the attitude control torque and G consists of the articu- 
lated joint torques. As in the previous cases, the attitude and angular velocity are 
measured at the attitude control torque actuator locations. Define 



a 


OJ 

Vp = 

e 

j Ur — 

e 


(48) 


In [7], the following control law was presented: 

_ + h{P + ^)Gpi + upis} 0 

[ 0 Gp2. 

where Gpi G Gp 2 G G jfj(^e+3)x(ne+3) symmetric positive defi- 

nite. It was shown in [7] that, if 0 < Am(Cpi) < there exists a unique equilibrium 
in the physical space, and that this equilibrium is globally asymptotically stable. 
Thus, the control law provides globally asymptotically stable large-angle maneuver- 
ing and posture control which can be used prior to capture to attain the desired 
end-effector positioning. 


Up GrPr (^ 9 ) 


4 Concluding Remarks 

Stabilization and control issues in capture and control of non-cooperative space ob- 
jects were addressed. The capture spacecraft was assumed to consist of a multi-link 
manipulator characterized by a base body with articulated appendages. The main 
problem of post-capture stabilization of the combined spacecraft/captured object 
was considered using proportional (variable-thrust) thrusters and combinations of 
thrusters and torque actuators. Control laws that asymptotically eliminate all rota- 
tional and translational velocities as well as all elastic motion, were presented. The 
control laws provide global asymptotic stability in spite of unknown mass-inertia 
properties of the captured object, as well as nonlinearities, unmodeled elastic modes, 
and uncertainties. The stability was shown to be preserved in the presence of (0, oo)- 
sector actuator nonlinearities (such as saturation). Combined attitude, position, and 
posture control was also considered using nonlinear quaternion-based feedback con- 
trol laws, which were shown to provide global asymptotic stability. An important 
issue that needs to be addressed is control during the capture stage when unknown 
contact forces are exerted on the capture spacecraft. Because the stabilizing control 
laws would be in effect during capture, the contact forces would act as bounded dis- 
turbances on the closed-loop system and would likely result in bounded rotational 
and translational velocities of the spacecraft, which would eventually reduce to zero 
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after capture is complete. However, this issue needs to be investigated further us- 
ing concepts such as input-to-state stability. It is also desirable to investigate the 
nse of on/off (rather than proportional) thrusters, which would be expected yield 
ultimate boundedness (e.g., a limit cycle) rather than asymptotic stability. Further 
research is also needed on control gain design and optimization, estimation methods 
for the object’s state, contact forces, and possibly mass-inertia properties, control 
methods for safe transportation of the object to another location, as well as detailed 
simulation studies and laboratory experiments. 
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